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The Standard Model of particle physics assumes that the so-called fundamen- 
tal constants are universal and unchanging. Absorption lines arising in molec- 
ular clouds along quasar sightlines offer a precise test for variations in the 
proton-to-electron mass ratio, ^^, over cosmological time and distance scales. 
The inversion transitions of ammonia are particularly sensitive to ^ com- 
pared to molecular rotational transitions. Comparing the available ammo- 
nia spectra observed towards the quasar B0218+357 with new, high-quality 
rotational spectra, we present the first detailed measurement of fi with this 
technique, limiting relative deviations from the laboratory value to | A/x//i| < 
1.8 X 10~^ (95% confidence level) at approximately half the Universe's current 
age - the strongest astrophysical constraint to date. Higher-quality ammonia 
observations will reduce both the statistical and systematic uncertainties in 
these measurements. 

The Standard Model of particle physics assumes that the fundamental constants of nature 
(or, at least, their low-energy limits) are the same everywhere and at every epoch in the Uni- 
verse. However, it cannot itself justify this assumption, nor can it predict their values. Our 
confidence in their constancy stems from Earth-bound laboratory experiments conducted over 
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human time-scales. Extrapolating to the entire Universe seems unwise, especially consider- 
ing that the physics driving the Universe's accelerating expansion, labeled 'Dark Energy', is 
completely unknown. Nevertheless, the Earth-bound experiments achieve impressive preci- 
sion: time-variations in the fine-structure constant, a = /he, which measures the strength 
of electromagnetism, are limited io a/a = (—1.6 ± 2.3) x 10^^'' yr^^ (i), while those in the 
proton-to-electron mass ratio, fi = m^/m^ - effectively the ratio of the strong and electro-weak 
scales - are limited to fi/ fi = (+1.6 ± 1.7) x 10~^^yr~^ (|2|). Still, more dramatic variations 
might have occurred over the 13-14 Gyr history of the Universe and the residual variations 
in our small space-time region might remain undetectably small. It is therefore imperative to 
measure the constants over cosmological time- and distance- scales. 

Variations in n and/or a would manifest themselves as shifts in the transition energies of 
atoms and molecules. By comparing transition energies registered in spectra of astronomical 
objects with laboratory values, possible variations can, in principle, be probed over our entire 
observable Universe and through most of its history. Due to the narrowness of the spectral 
features involved, absorption lines arising in gas clouds along lines of sight to background 
quasars are currently our most precise cosmological probes. For example, by comparing various 
heavy element electromagnetic resonance transitions in optical quasar spectra, 5-a evidence has 
emerged for variations in a of ~5 parts in 10^ at redshifts 0.2 < z < 4.2 (|J1|?1I2]I3). Although 
a more recent statistical sample found no variation (7), errors in the analysis prevent reliable 
interpretation of those results (019!), leaving open the possibility of a varying a. 

Similarly, tentative (3-cr) evidence for a fractional variation in fi of ~20 x 10"^ has recently 
come from 2 quasar spectra containing many ultraviolet (UV) H2 transitions at redshift z ~ 
2.8 (17011771) . Comparison of UV heavy-element resonance lines and Hi 21-cm absorption is 
sensitive to a^/x and, assuming a to be constant, has yielded indirect null constraints on fi- 
variation, albeit with slightly worse precision than the direct H2 method (72 ). 

An alternative method for measuring /i at high redshift, suggested recently by Flambaum & 
Kozlov (173|) . is to utilize the sensitivity to yU- variation of the ammonia inversion transitions (fT?]) 
near 24 GHz. A shift in their frequencies due to a varying ^ can be discerned from the cosmo- 
logical redshift by comparing them to transitions with lower sensitivity to ji. Good candidates 
for comparison are rotational transitions of molecules such as CO, HCO+ and HCN because (i) 
their transition frequencies depend mainly on /i and not on other fundamental quantities (e.g. a); 
(ii) they are simple molecules, commonly detected in the interstellar medium; (iii) their rest fre- 
quencies (80-200 GHz) are not vastly dissimilar to the NH3 transitions' (cf. comparison of UV 
and 21-cm), thus reducing possible effects due to frequency-dependent spatial structure in the 
background quasar's emission. This final point is discussed in more detail below. 

For rotational and NH3 inversion transitions we may write the apparent change in velocity 
or redshift of an absorption line due to a variation in fi as 

c I + z n 

where Ki is the sensitivity of transition i to /i and /S^ji/ fi = (/x^ — jj^iab) / f^iab for /iiab and 
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the values of /i in the laboratory and absorption cloud at redshift z, respectively. All rotational 
transitions have Ki = 1, so comparing them with each other provides no constraint on A/i//i. 
However, the sensitivity of the NH3 inversion transitions is strongly enhanced, Ki ^ 4.46 (|73l) . 
That is, as ji varies the NH3 transitions shift relative to the rotational transitions. 

Only one quasar absorption system displaying NH3 absorption is currently known, that 
at 2; = 0.68466 towards quasar B0218+357. From redshift uncertainty estimates for NH3, 
CO, HCO+ and HCN in the literature, (fTj)) crudely estimated the precision achievable to be 
5(A/i//i) ^ 2 X 10~^. However, (T3) cautioned that a proper measurement of A^/ fi from NH3 
would require detailed, simultaneous fits to all the molecular transitions and that significantly 
better precision may be possible. Here we make the first detailed measurement of ^ using the 
NH3 inversion transitions by comparison with HCO+ and HCN rotational transitions. 

The only published NH3 absorption spectra are those for the (J,K) = (1,1), (2,2) and (3,3) 
inversion transitions reported by (i75l) towards B02 18+357, reproduced in Fig. 1. The channel 
spacing is 1.67kms^^ for(l,l) & (2,2) and3.3kms~^ for (3,3). The spectra are normalized by 
a low-order fit to the quasar continuum. See (|75l) for observational and data reduction details. 
The signal-to-noise ratio (SNR) for the flux is very high, ~ 1000 per channel, but since <1% of 
the continuum is absorbed, the effective SNR for the optical depth is quite low. 

HCO"*" and HCN(l-2) absorption towards B0218+357 was discovered more than a decade 
ago (16^ . New, high resolution (^0.9 km s^^ per channel), high SNR (~100 per channel) obser- 
vations of these lines were recently undertaken with the Plateau de Bure Interferometer, France; 
see (|77|) . Fig. 1 shows both spectra normalized by fits to the quasar continuum. 

All spectra were registered to the heliocentric reference frame; possible errors in this proce- 
dure are discussed in {17) and shown to be negligible. 

Spectra representing the 1-cr uncertainty in normalized flux per channel were constructed 
for all the molecular spectra by calculating the root-mean- square (RMS) flux variations in the 
continuum portions of each transition. Since no large differences were observed either side of 
the absorption for any transition, a simple constant error model was adopted. 

As equation 1 states, the signature of a varying ji would be a velocity shift between the rota- 
tional and NH3 inversion transitions. Complicating the measurement of any shift is the 'velocity 
structure' evident in Fig. 1: the profiles comprise absorption from many gas clouds, all asso- 
ciated with the absorbing galaxy but nevertheless moving at different velocities. The number 
and velocity distribution of these 'velocity components' are unknown and must be determined 
from the data themselves. Each fitted velocity component is represented by a Gaussian profile 
parametrized by its optical depth, Doppler width and redshift. The best-fitting parameter val- 
ues are determined with a -minimization code, VPFIT ^T8\), designed specifically for fitting 
quasar absorption lines. To determine the statistically preferred velocity structure, the best fit 
(i.e. minimized) values of per degree of freedom, xl, are compared for several fits with dif- 
ferent velocity structures. That with the lowest xl is taken as the fiducial one (i.e. an F-test to 
discriminate between models). 

Measuring Ayu/yU requires the assumption that the velocity structure is the same in all transi- 
tions. This does not mean that the ratio of optical depths of corresponding velocity components 
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in different transitions must be constant across the profile. Rather, it means that the number 
and velocity distribution of components are assumed to be the same in different transitions. We 
discuss this assumption below but, since it must be made eventually, in practice the velocity 
structure was determined by tieing together the redshifts of corresponding velocity components 
in different transitions. The high SNR rotational spectra clearly place the strongest constraints 
on the velocity structure but the NH3 spectra must be included to measure A/i//i. The Doppler 
widths of corresponding components were also tied together, effectively assuming a turbulent 
broadening mechanism. 

Figure 1 shows the fiducial 8-component fit. Note that the detailed hyperfine structure of 
the HCN(l-2) and NH3 transitions is reflected in each velocity component. The relative hy- 
perfine level populations were fixed by assuming local thermodynamic equilibrium (LTE). The 
laboratory data used in the fits are tabulated in (I77|) . Given this fit, determining A/i//i is straight- 
forward: a single additional free parameter is introduced for the entire absorption system which 
shifts all the velocity components of each transition according to its /^-coefficient (equation 1). 
All parameters in the fit, including the single value of A/i/yU, are varied by VPFIT to minimize 
X^. The best-fit value is A/i/yU = (+0.74 ± 0.47) x 10^^, corresponding to a (statistically in- 
significant) velocity shift between the NH3 and rotational transitions of 0.77 ±0.49 km s^^. The 
l-cr error quoted here - formed from the diagonal terms of the final parameter co variance ma- 
trix - derives entirely from the photon statistics of the absorption spectra. A different (though 
intimately related) approach to determining A/x/yU and its error is discussed in (|77l) : it provides 
the same result. 

Fitting too few velocity components causes large systematic errors in such analyses dP])- 
A reliable measurement of Afi/ fi can only be derived from fits replicating all the statistically 
significant structure in the absorption profiles. Therefore, the fiducial velocity structure must be 
the statistically preferred one. When fitting many transitions simultaneously, this may be more 
complicated than the one preferred by human eye. The simplest objective method to achieve 
this is demonstrated in Fig. 2 which shows the decrease in as increasingly complex velocity 
structures are fitted. When the fit is too simple to adequately describe the data, quite different 
values of A/i/zi are found. On the other hand, the 9-component 'over-fitted' case provides a 
value and error very similar to the fiducial 8-component model. 

The consistency of the velocity structures in the two highest SNR transitions, HCO+ and 
HCN(l-2), was tested by fitting those transitions independently. Again, different fits with dif- 
ferent velocity structures were compared to determine the statistically preferred one. Both tran- 
sitions are best fit by velocity structures similar to that in Fig. 1 [see (i7)], providing some 
confidence that they can meaningfully be fitted simultaneously. These independent velocity 
structures were applied to the NH3 transitions and new values of Ayu//i derived. When fit- 
ting only HC0+(l-2) and NH3, A/i//i = (+0.67 ± 0.51) x 10^^ for HCN(l-2) and NH3, 
A/i//i = (+0.88 ± 0.51) X 10^^. Neither value substantially deviates from our fiducial one. 
Note that the marginal increase in the l-cr error when using a single rotational transition indi- 
cates that the NH3 spectra limit the statistical uncertainty. 

To consider potential systematic uncertainties, it is important to recall our main assumption: 
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the velocity components constituting the absorption profiles have the same redshifts in different 
transitions. While the HCO+ and HCN(l-2) velocity structures are evidently similar enough 
for measuring Ayu//i, the NH3 spectra have too low SNR for a direct comparison. And since the 
observed frequencies of the NH3 (~ 14 GHz) and rotational (~ 106 GHz) transitions are some- 
what different, it is possible that, if the background source morphology is frequency dependent, 
some NH3 components might arise along slightly different sight-lines to those components in 
the rotational profiles. This is considered presently. 

B0218+357 is a z = 0.944 BL Lac object {19) lensed by a nearly face-on Sa/Sab z = 0.68 
galaxy (20) in which the absorption occurs. Two lensed images, A and B, separated by 334 mas, 
straddle the lensing galaxy's center with image B much closer to the center. An Einstein ring 
with diameter ~300mas, centered near image B, has also been identified (|27T) . B0218+357 
itself has a core-jet morphology with an unresolved [<1 x 1 mas or <7x 7pc (22)] flat-spectrum 
core dominating the observed 8.4 GHz emission (l23i) . The jet has a knotty structure extending 
over ~10 X 10 mas and, like other jets, is expected to have a steep spectrum. 

Various absorption lines have been detected in the z = 0.68466 absorber, from Hl21-cm 
and OH below (rest-frame) 2 GHz (|2l|25l), through six H2CO transitions at 5-150 GHz (|2gl) . 
to H2O at 557 GHz (|27l) to name but a few. Furthermore, the molecular absorption arises only 
towards image A (|23l29ll30ll23l) . The flat-spectrum core should completely dominate at high 
frequencies; that H2CO and H2O absorb most of the total high-frequency continuum therefore 
implies that at least those transitions only arise towards the core. All the observed molecu- 
lar transitions have consistent velocity structures (though most spectra have poorer resolution 
and/or SNR than those studied here). Thus, the most important velocity components in all 
transitions evidently arise towards image A's flat-spectrum, compact core (|75ll30|) . 

Nevertheless, indirect evidence suggests that the molecular clouds do not completely cover 
the background source (137] [771 ). If the covering fraction is different for the rotational and NH3 
transitions, some velocity components may appear in one and not the other. Similar problems 
may arise because some HCO+ and HCN(l-2) velocity components may be optically thick (fTTl) . 
The spurious shifts in A|u//i these effects may cause are difficult to estimate in general but 
in (17) we conduct several fits in which different combinations of NH3 velocity components 
are removed, providing an estimate of ±0.7 x 10"^. Another potential systematic error is our 
assumption of LTE for the HCN(l-2) and NH3 hyperfine structure populations. As mentioned 
above, removing HCN(l-2) from the analysis barely changes the measured Afi/ fi. Removing 
different parts of the hyperfine structure from the NH3 transitions results in maximum deviations 
of ±0.3 X 10"^ from our fiducial value of Afi/fi (17). 

Combining these two potential systematic errors in quadrature, we obtain A/i//i = (±0.74± 
0.47stat ± 0.76sys) X 10^^, providing no evidence for cosmological variations in /i. The NH3 
spectra currently set both the statistical and systematic errors. While, clearly, the SNR and 
resolution directly determine the former, they also indirectly influence the latter: the velocity 
structure of higher quality NH3 spectra could be more directly compared with the rotational 
profiles and limits on non-LTE hyperfine structure anomalies could be constrained by the data 
themselves. That is, with improved NH3 spectra, both the statistical and systematic error com- 
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ponents can be improved. Nevertheless, until the NH3 data are improved, our final result is a 
2-cr limit on variation in fi from this single absorber: |Ayu/yu| < 1.8 x 10^^. This corresponds 
to < 1.9 km s^^ shift between the NH3 and rotational transitions. 

Our new value of A/i//i seems inconsistent with the current tentative evidence for /i-variation 
from H2 absorption at 2; ~ 2.8, A/i//i = (+24.4 ± 5.9) x 10~^ (|72]). However, reliable com- 
parison is difficult because cosmological time- and/or space-variations in /i remain poorly con- 
strained. Clearly, a statistical sample from both techniques, covering a wide redshift range and 
with detailed assessment of systematic effects, is highly desirable. 

Assuming /x varies linearly with time, our measurement corresponds to a drift of = 
(—1.2 ± 0.8stat ± l-2sys) X 10~^^ yr^^- However, this assumption is only a convenient means for 
comparison with current limits from laboratory atomic clocks, /i//i= (+1.6±1.7)xl0^^^yr^^ 
(I2I); it is not motivated by any physical (necessarily untested) varying-/i theory. 

With so few measurements of jj, distributed throughout the Universe, each new measurement 
is an invaluable test of the most basic - and theoretically unjustifiable - assumptions in the Stan- 
dard Model: that the laws of physics are universal and unchanging. The precision demonstrated 
here highlights the importance of discovering many more molecular absorbers to our knowledge 
of fundamental physics. 
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Figure 1: Spectra of the molecular transitions used in this study registered to a heliocentric 
velocity scale centered on z = 0.68466. The nominal observed frequencies are noted in each 
panel. The data, normalized by fits to their continua, are plotted as black histograms. Tick- 
marks above the spectra show the positions of velocity components in our fiducial 8-component 
fit (solid line following the data). Note that the HCN and NH3 transitions have complex hy- 
perfine structure reflected in each velocity component; the tick marks show the position of the 
strongest hyperfine component in LTE [see ((77l) 1. Residuals between the fit and data, normalized 
by the (constant) error array, are plotted above the spectra, bracketed by horizontal lines repre- 
senting the ±1 cr level. The fit contains 57 free parameters: an optical depth for each component 
in each transition (5x8 parameters) plus a Doppler width and redshift for each component (8+8 
parameters) and a single value of A/x/ /i. The fitted line parameters are tabulated in (|77|) . 
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Figure 2: Variation in A/i/yU and per degree of freedom, xl^ of different velocity structures 
characterized by the number of fitted absorption components, xl is defined as x^/i^ = Jlf"^ [dj — 
m{j)]/aj for dj the j^^ data value with variance (t| and model value m{j). The sum is over 
all A^'d = 223 data points; u = — A^par for A^par free model parameters. Our fiducial 8- 
component (A^par = 57) result is highlighted with square points. Different components were 
added/removed to/from the fiducial fit to form each initial velocity structure and VPFIT was run 
again to minimize by varying all free parameters. Two different initial fits with 6 components 
and three fits with 7 components were possible; the different results are offset in the plot for 
clarity in these cases. Large xl values for < 6 components indicate that those fits are not 
statistically acceptable. Of the remaining fits, the 8-component fit has the lowest xl- Note 
that the 9-component fit has a smaller x'^ (because more parameters are being fitted) but a 
marginally higher xl, indicating that it is less statistically preferred than the 8-component fit. 
Only statistical error bars on Ayu//i are shown; see text for discussion about systematic errors. 
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SUPPORTING ONLINE MATERIAL 



Materials and methods 

New observations of HCO+ and HCN(l-2) 

The HCO+ and HCN(l-2) spectra were obtained with the Plateau de Bure Interferometer 
(PdBI) with the 3-mm receivers, the HC0+(l-2) line being redshifted to 105.882 GHz and 
the HCN(l-2) line to 105.221 GHz. Five or six antennas were used in compact configurations, 
and the two images, A and B, of B02 1 8+357 were not resolved from each other. The continuum 
flux density of B02 18+357 was estimated to be ~0.5 Jy. 

The spectra presented here are averages of spectra taken on several observing runs: 28 and 
29 June 2005, 15 July 2006 and 23 March 2007 for the HCO+ data; 5 July 2005, 19 August 
2006 and 23 March 2007 for the HCN data. The integration time varied from 30-min to ~2- 
hr. Within the SNR of individual integrations, we do not see time variation in the absorption 
profiles between 2005 and 2007. Note that in March 2007 the HCO+ and HCN data were 
obtained simultaneously with dual polarization receivers. 

The frequency dependence of the instrument across the spectral bandwidth (80 MHz with 
256 channels of width ~0.3 MHz) was calibrated by observing a strong radio quasar for several 
minutes. Its amplitude and phase as a function of radio frequency were solved for by fitting 
high order polynomials (10 and 20, respectively). The RMS residuals of the bandpass cali- 
bration were typically less than 1-2%. The visibilities were self-calibrated with respect to the 
continuum emission of B02 18+357 itself: the phase was referenced to the barycenter of the 
continuum emission and the amplitude was normalized to the total continuum flux. The spectra 
observed with the different baselines of the interferometer and at different epochs were then co- 
added. Third-degree baselines (derived on channels free of absorption) were subtracted from 
these final spectra to remove small residuals in the radio bandpass calibration, instrument drifts 
and/or atmospheric fluctuations. 

Laboratory molecular data 

Table S 1 provides all the laboratory molecular data used in our analysis. Laboratory frequencies 
for HCO+ and HCN(l-2) were taken from (SI) and (52) respectively, as recommended by 
(1531) . Frequencies for the NH3 inversion lines were taken from (15?)) . Including the hyperfine 
structure of the HCN(l-2) and NH3 inversion transitions when fitting the absorption spectra is 
very important. The relative optical depths of the HCN(l-2) hyperfine components, assumed to 
be in LTE, were taken from (S3 ) while those for the NH3 lines were computed using the formula 
presented by (S4!i . 
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Alternative minimization scheme 

The simplest, fastest and most robust method of determining /S.^/ fi from a multi-component, 
simultaneous fit to different transitions is simply to include A/i//i as an additional free param- 
eter in the minimization process. This was our approach in the main text. An alternative 
approach is to conduct many simultaneous fits to the absorption spectra, each with a fixed input 
value of A/i//i but allowing all other free parameters to vary, to build up a 'y^ curve' - a plot 
of versus A/x/ /x. The best-fitting value of Afi/ fi is the one which produces the lowest value 
of Xmin- The functional form of implies that the curve should be near parabolic and 
smooth on scales similar to the l-a error on A/x/yU. The error can therefore be determined from 
the width of the curve near the minimum; the curve's smoothness is obviously essential to 
this procedure - see cf. (J}. 

Figure SI shows the y^ curve for our fiducial 8-component fit to the rotational and NH3 
inversion spectra. For convenience, we plot Ax^ = y^ — y^^^ normalized by the lowest value 
of yl - i.e. y^ per degree of freedom - in the curve, X^min- This representation allows the 
\-a error in Ajj,/ jj, to be determined directly from the width of the curve at Ax^/x^min = 1. 
The result from this process is A/x/yU = (+0.75 ± 0.45) x 10^^, in very close agreement 
with the value determined simply by fitting A/i//i as an additional free parameter, A/i//i = 
(+0.74 ±0.47) X 10"^ 

Errors from bandpass calibration and baseline stability? 

Possible observational systematic uncertainties include bandpass calibration errors and vari- 
able spectral baselines. The former may affect the PdBI (i.e. iterferometric) rotational spectra 
while the latter will be most severe for the single-dish NH3 spectra which were obtain with 
the Effelsberg 100-m Ratio Telescope in Germany. Both these effects could potentially lead 
to small ripples in the absorption spectra. As mentioned above, the RMS bandpass calibration 
errors for the new rotational spectra are <2%. Additionally, the fact that A/i//i is insensitive to 
the rotational transition used in the analysis suggests that any residual ripples in the HCO^ or 
HCN(l-2) spectra are unimportant. The poorer SNR (on the optical depth) of the NH3 transi- 
tions makes individually removing them a less stringent test. Nevertheless, removing NH3 (1,1) 
yields A^u/yU = (+0.67 ± 0.62), removing NH3 (2,2) gives (+0.65 ± 0.62) and eliminating 
NH3(3,3) provides (+0.84 ± 0.49), all in units of 10"*^. Wald-Wolfowitz runs tests on the 
unabsorbed portions of the NH3 spectra also reveal no evidence for significant ripples. There- 
fore, there is no evidence for significant systematic effects due to bandpass calibration and/or 
baseline stability problems. 

Alternative velocity structures 

As described in the main text, we fitted several different velocity structures in our search for 
the one with the lowest y^, i.e. the fiducial 8-component model plotted in Fig. 1. For each 
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velocity structure, VPFIT was run again to minimize by varying all free parameters of the fit. 
These alternative fits provided the black points in Fig. 2. For comparison, Fig. S2 plots these 
alternative structures in the same way as Fig. 1 . By inspecting the residual spectra (plotted above 
each transition), it is easy to see why the fits with fewer than 8 components result in larger values 
of xl compared to the fiducial model. The residuals take many -pixel excursions outside the l-cr 
ranges at similar velocities in different transitions - a sure sign that real structure exists in the 
absorption profiles which is left unaccounted for by the fit. For example, the human eye readily 
picks out 5 or 6 components in the rotational profiles but the two different 6-component models 
produce very poor fits (two right-hand panels in middle row) with xl > 1-6 in Fig. 2. On the 
other hand, the 9-component model provides "too good" a fit, as demonstrated by its marginally 
higher xl value in Fig. 2 compared to the fiducial model. The best-fit line parameters for the 
fiducial 8-component model are given in Table S2. 

The high SNR of the HCO+ and HCN(l-2) spectra implies that our determination of the 
velocity structure is driven by those transitions and not the NH3 inversion lines. A good cross- 
check on the central assumption in our analysis - that the velocity structure is the same in differ- 
ent transitions - can therefore be made by fitting the HCO^ and HCN(l-2) spectra separately. 
Fig. S3 compares our independent fits of these rotational lines. Only the velocity structures 
providing the lowest xl values are plotted. Seven components provide the lowest-x^ fit for the 
HCO+ data while 8 components are required to fit the HCN spectrum. However, as for our fidu- 
cial simultaneous fit in Fig. 1, the bluest component in HCN is not statistically required by the 
HCO+ data. That is, both the HCN and HCO+ data can be individually fit by the same number 
of components as used in the simultaneous fit. Also note the general agreement between the 
velocities at which these components are best fitted. In this sense the velocity structures of the 
two high SNR rotational transitions are consistent. 

Errors from morphology and variability of B0218+357? 

Quasars are not point sources at radio wavelengths but instead show sometimes complex mor- 
phological structure. Indeed, each lensed image of B02 18+357 shows a core-jet morphol- 
ogy. If the absorbing molecular clouds at z = 0.68466 have smaller angular extent than the 
background radio source, and the source's morphology varies between 24 and 178 GHz (in the 
absorber's reference frame), then the NH3 and rotational transitions may arise along slightly 
different sight-lines. This may, in principle, lead to systematic errors in the determination of 
A/i/ fi. 

As discussed in the main text, all the molecular absorption lines studied here are very likely 
to arise only towards the core of image A. Its angular extent, <1 x 1 mas, places an upper limit 
on the spatial extent of the absorbing region of just <7 x 7pc (175] 13011221) . This is a typical 
size for Galactic molecular clouds (S5). Nevertheless, if the core is close to the upper limit in 
angular extent and/or if the absorbing clouds are smaller than ~5 pc, then it is possible that the 
NH3 and rotational profiles may have different velocity structures. The similarity of the velocity 
structures of a variety of transitions - covering a larger frequency range than just those studied 
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here - suggests this may be of little concern. 

However, evidence for small-scale structure of the absorbing clouds has been found in (|37T) . 
By examining the optical obscuration and reddening of image A by dust expected to reside in 
the absorbing clouds, (31) established that not all of image A's optical emission is obscured 
and that its unobscured portion is not even substantially reddened. Therefore, assuming that 
the main component of optical emission originates from image A's core (just as for the radio 
emission), the absorbing clouds must not completely cover the background emission. Since 
the emission region's size should vary with frequency, the NH3 profiles may contain a different 
number of absorbing velocity components compared to the rotational profiles (though some 
components may still be in common). 

The flux variability of B02 18+357 may also be of potential concern. Again, if the absorbing 
clouds do not cover the core of image A and its morphology changes with time, then observa- 
tions at different epochs may probe slightly different sight-lines. As detailed above, the HCO+ 
and HCN(l-2) data were obtained over several runs between June 2005 and March 2007 with 
HCO"*" usually observed up to 1 month before HCN. The NH3 data were collected in August 
2001 and June 2002. 

It is difficult to use the available data to test specifically for each of these effects. Neverthe- 
less, the fact that the low SNR NH3 profiles at least resemble the higher SNR rotational ones - 
particularly the velocity of the peak optical depth - allows a general quantitative assessment of 
the likely effect on the Afi/ fi measurement. One test is to remove specific (groups of) velocity 
components from the fits to the NH3 spectra while keeping the fiducial 8-component model for 
the rotational lines. Upon removing a single velocity component - that is, fitting 7 components 
to the NH3 spectra - the different possible fits yield an RMS deviation of 0.5 x 10^^ from our 
fiducial A/i/yU. This is dominated by the results of removing either component 3 or component 
6 (see Fig. 1), which give deviations of ~ 0.7 x 10~^. Removing selected pairs of components 
(e.g. 1 & 3 or 5 & 7 etc.) or selected combinations of three components yield similar RMS 
deviations of ^ 0.5 x 10"^; these are again dominated by cases where either components 3 or 
6 are removed. Another test is to fit the NH3 profiles with as few components as is statistically 
acceptable (i.e. ~ 1) - in our case just two, components 3 & 6 - while keeping the full 
velocity structure in the rotational lines. For this test we relaxed the constraint that the widths 
of the rotational and NH3 components should be the same. We find a deviation of 0.45 x 10^^ 
from the fiducial value to A^/^u = (+0.29 ± 0.39) x 10"^ 

Given these results, we adopt a representative systematic error on Ayu//i of ±0.7 x 10~^ 
due to possible differences in the velocity structure of the rotational and NH3 spectra. Note 
that if higher resolution, higher SNR NH3 spectra were obtained, a more direct comparison of 
the NH3 velocity structure with that of the rotational transitions could be made, similar to the 
comparison of the two rotational transitions made above. Therefore, high quality NH3 spectra 
should reduce the systematic error calculated here. 
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Optical depths 



Neither the PdBI rotational observations nor the Effelsberg NH3 observations resolved the two 
lensed images from each other. Thus, even though the molecular absorption is very likely to 
arise only towards the core of image A, continuum flux from image B, the jet of image A 
and possibly the Einstein ring will also contribute to the spectra studied here. The flat-spectrum 
lensed cores of images A and B are thought to dominate the continuum flux beyond 14 GHz and 
image A contributes >70% of the flux ([S6][75l l?D|). Therefore, the fraction of image A's core 
continuum absorbed by the NH3 and rotational lines is about a factor of 1 .4 larger than shown 
in Fig. 1. Thus, if some velocity components truly absorb all of image A's core continuum (at 
the observed spectral resolution), they will appear flat-bottomed but will not extend down to the 
zero flux level. 

For simplicity, we did not attempt to take this potential effect into account in the profile 
fits in Fig. 1 . It is clearly not a problem for the NH3 inversion lines since they only absorb 
~1% of the total continuum flux. However, it is possible that some velocity components in 
HCO"*" and HCN(l-2) are near saturation (see below), even though we do not observe flat- 
bottomed profiles. We tested the effect of this scenario on A/i/|U in the following way: 30% of 
the flux was subtracted from the rotational spectra and the continuum levels were appropriately 
renormalized. Instead of Gaussian absorption components we fitted Voigt profiles to allow 
saturation effects in the fit. We conducted several fits with a range of quantum mechanical 
damping coefficients to characterize the onset of damping wings in the Voigt profiles. All fits 
gave results negligibly different to our fiducial result, indicating that, under the assumption that 
our fitted velocity structure reflects the real velocity distribution of absorbing clouds, saturation 
effects do not affect our measurement of Ayu//i. 

It is likely that some of the stronger components in the HCO+ and HCN(l-2) profiles are 
optically thick. Precise optical depths are difficult to determine because the fraction of the 
background source covered by the absorption clouds is also unknown. Wiklind & Combes (1^ 
found optical depths of > 1.9 for HCN(l-2) and > 1.1 for HC0+(l-2). A similar conclusion 
was reached for HC0+(l-2) by (|30l) . If several velocity components with very similar velocities 
(Av<2kms~^) and different optical depths exist in the rotational spectra but are fitted as a 
single component, then some spurious shifts with respect to the optically thin NH3 profile might 
be measured. The potential systematic error on A|u//i is therefore the same as that caused by 
uncertainty in the NH3 velocity structure, as quantified above. 

Errors from hyperfine structure? 

All the fits discussed in the main text and above assumed that the relative optical depths of the 
HCN(l-2) and NH3 hyperfine structure components were those expected in LTE - see Table 
SI. However, non-LTE hyperfine structure effects have been observed in HCN(O-l) emission 
lines from Galactic interstellar clouds (|57l) . Nevertheless, if LTE conditions do not prevail for 
the HCN(l-2) velocity components, then the effect on Afi/ fi seems to be small since remov- 
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ing HCN(l-2) from the analysis does not change our measured A^/ ^ substantially (see main 
text). However, removing the NH3 transitions is not possible if we are still to measure Ayu/yU, 
nor are the NH3 spectra of high enough quality to directly determine the hyperfine structure 
populations; the potential effect of non-LTE hyperfine structure must be assessed differently. 

NH3 hyperfine structure 'anomalies' have previously been found in Galactic sources both 
in absorption ([S^ and emission ([S9l) . The non-LTE conditions cause the AFi 7^ 'satel- 
lite' components with a given quadrupole quantum number, Fi, to have different optical depth 
depending on the change in Fi during the transition. Since complementary (i.e. Fi=l-2 and 
Fi=2-1) transitions straddle the main (AFi = 0) hyperfine components, non-LTE conditions 
cause asymmetries in the optical depth profile of any given velocity component. A simple test 
to gauge the maximum effect this may have on A/x/yU is to remove either the high or low fre- 
quency satellite lines from the NH3 fits. This causes RMS deviations from our fiducial Afi/fx 
of 0.3 X 10~^. We adopt this as a second component to our systematic error budget in the main 
text. Note that with higher resolution, higher SNR NH3 spectra it should be possible to con- 
duct consistency checks for strong non-LTE hyperfine structure effects, thereby reducing the 
possible systematic error. 



Reference frames 

All the spectra studied here were shifted to the heliocentric rest frame so that any velocity shifts 
between them could be attributed to a varying /i. 

Depending on the format in which raw radio telescope data are recorded, passed to the 
data reduction software and processed, one or more small velocity shifts may be applied to 
them. These typically convert between the local standard of rest (LSR) and heliocentric frames. 
However, the Doppler correction for the radial velocity, 

2 I 2 ' ^ 
C '^obs + ^0 

is approximated in many data reduction software packages (for traditional reasons, not ones of 
simplicity or processing speed etc.) by 

~ (3) 

which applies when ^rad/c 1. Here uq is the emitted frequency and z/obs is the observed 
frequency. Since we are interested in small shifts between different transitions with a precision 
of 6{Afj,/ fi) ~ 0.1 X 10^^, corresponding to a velocity shift of just 0.1 kms^\ we must check 
whether the above approximation is accurate enough so as not to cause systematic errors in 
A/i/ /i. The difference between the above two expressions for the Doppler correction is approx- 
imately (frad/c)^ for frad/c ^ 1. So, for a maximal heliocentric correction of ~30kms~\ 
the systematic error introduced by using the approximation is ~3ms~^. This corresponds to 
5(A/i//i) ~ 3 X 10^^, well below our current statistical precision. 
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Figure SI: curve for the fiducial 8-component model fit. Note its (required) smoothness 
and (expected) near-parabolic shape. The right-hand vertical scale shows xl which takes a 
minimum value of xlmm — 1-07125. The left-hand vertical scale shows Ax'^/xlmin^ thereby 
allowing the l-a error in A/x/yU to be immediately read off the graph (dotted black lines). In 
practice, we measured the position at which takes its minimum and the curve's width by 
fitting a parabola to the central points; the fit, marked by the black solid line, demonstrates 
how close to parabolic that part of the curve is. The black circle and error-bar represent the 
result from this fit, Afi/ fi = (+0.75 ± 0.45) x 10^®, which closely matches our fiducial value 
represented by the red/grey square and error-bar. 
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Component 9 added 



NHgCl, 1) 



Velocity around z = 0. 66466 [km s"'] 
Component 8 removed 



HC0-(l-2) 



Lr-! NHgCa.a) 



NH3(3,3) 

Velocity around z = 0. 68466 [km s"'] 
Components 3 4&B removed 



HC0*(l-2) 



HCN (1-2) 



Velocity around z = 0. 68466 [km s '] 




Component 3 removed 



HCN (1-2) 



NHgCl, 1) 



Velocity around z=0. 68466 [km s"' 
Components 4 &: 8 removed 



NHg (3, 3) 

Velocity around z=0. 68466 [km s~'] 
Components 3, 4, 7 & 8 removed 



HCN (1-2) 




Component 4 removed 



^^.J' NH3(1,1) 



NH3(3,3) 



Velocity around z = 0. 68466 [km s"'] 
Components 3 & 4 removed 



HCO* (1-2) 



HCN (1-3) 



~\-J~ NHg (2, 2) 



NHg (3, 3) 

Velocity around z = 0. 68466 [km s"'] 



Velocity around z=0. 68466 [km s 



Figure S2: The different velocity structures attempted in search of the statistically preferred 
one. The values of A/i/yU and xl corresponding to each fit are represented by the black points 
in Fig. 2. Components were removed/added from/to the fiducial fit to form each initial velocity 
structure and vpfit was run again to minimize x^- For the 9-component model, the additional 
component was added redwards of all other components. The residual spectra show how poor 
the fits are with < 6 components. 
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Figure S3: Independent fits to the two rotational transitions studied here. Layout similar to 
Fig. 1. The data, tick-marks and residuals for HCN(l-2), marked in lighter grey lines, are 
offset below those for HCO+ for clarity. The formal l-a statistical error in the redshift of each 
component (converted to velocity) is represented by the horizontal bars across the tick marks. 
The error bar for each component is plotted higher than the one to its left for clarity. There is 
broad agreement between the positions of the velocity components. The bluest component in 
HCN(l-2) is not statistically required by the HC0+(l-2) data alone, similar to the simultaneous 
fit in Fig. 1. 
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Table SI: Laboratory frequencies of transitions used in our analysis. The second and third 
columns provide, respectively, the relevant quantum numbers for the lower and upper states in 
the following order: J, F, Fi. For HCO+(l-2) only J is relevant while Fi is not relevant for 
HCN(l-2). For the NH3 inversion transitions, J = K, F is the magnetic hyperfine quantum 
number and Fi is the quadrupole quantum number. The final column shows the relative optical 
depth of different hyperfine components, S, in LTE. E is normalized to unity for the strongest 
hyperfine component. 



Molecule 




Transition 


Frequency 


logs 




Lower 


Upper 


[GHz] 




HCO+ 




1 




2 


178.375065(50) 


0.0000 


HCN 




1,2 




2,2 


177.2596770(20) 


-0.7482 


HCN 




1,0 




2, 1 


177.2599230(20) 


-0.6233 


HCN 




1,1 




2,2 


177.2611100(20) 


-0.2711 


HCN 




1,2 




2,3 


177.2612230(20) 


0.0000 


HCN 




1,2 




2, 1 


177.2620122(11) 


-1.9243 


HCN 




1, 1 




2, 1 


177.2634450(20) 


-0.7482 


NH3 


1, 


1/2, 1 


1, 


1/2,0 


23.692926829(56) 


-0.7993 


NH3 


1, 


3/2, 1 


1, 


1/2,0 


23.692968829(38) 


-0.4983 


NHa 


1, 


1/2, 1 


1, 


3/2,2 


23.693872152(57) 


-0.7024 


NHa 


1, 


3/2, 1 


1, 


5/2,2 


23.693905112(58) 


-0.4472 


NH3 


1, 


3/2, 1 


1, 


3/2,2 


23.693914466(39) 


-1.4014 


NHa 


1, 


1/2, 1 


1, 


1/2, 1 


23.694459098(36) 


-1.1004 


NH3 


1, 


1/2, 1 


1, 


3/2, 1 


23.694470034(28) 


-1.4014 


NH3 


1, 


3/2,2 


1, 


5/2,2 


23.694470904(63) 


-1.1461 


NH3 


1, 


3/2,2 


1, 


3/2,2 


23.694480291(44) 


-0.1919 


NH3 


1, 


3/2, 1 


1, 


1/2, 1 


23.694501428(59) 


-1.4014 


NH3 


1, 


5/2,2 


1, 


5/2,2 


23.694505950(37) 


0.0000 


NH3 


1, 


3/2, 1 


1, 


3/2, 1 


23.694512322(37) 


-0.7024 


NH3 


1, 


5/2,2 


1, 


3/2,2 


23.694515319(47) 


-1.1461 


NH3 


1, 


3/2,2 


1, 


1/2, 1 


23.695067195(46) 


-0.7024 


NH3 


1, 


3/2,2 


1, 


3/2, 1 


23.695078206(59) 


-1.4014 


NH3 


1, 


5/2,2 


1, 


3/2, 1 


23.695113176(50) 


-0.4472 


NH3 


1, 


1/2,0 


1, 


1/2, 1 


23.696029719(48) 


-0.7993 


NH3 


1, 


1/2,0 


1, 


3/2, 1 


23.696040646(38) 


-0.4983 


NH3 


2, 


3/2,2 


2, 


3/2, 1 


23.720534302(33) 


-1.7058 


NH3 


2, 


5/2,2 


2, 


3/2, 1 


23.720575068(39) 


-0.7516 


NH3 


2, 


3/2,2 


2, 


1/2, 1 


23.720579876(34) 


-1.0068 


NH3 


2, 


5/2,2 


2, 


7/2,3 


23.721336248(58) 


-0.7570 


NH3 


2, 


3/2,2 


2, 


5/2,3 


23.721337256(57) 


-0.9119 


NH3 


2, 


5/2,2 


2, 


5/2,3 


23.721377951(106) 


-2.0580 


NH3 


2, 


1/2, 1 


2, 


3/2, 1 


23.722588837(70) 


-1.0068 
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Table S 1 : continued 



Molecule Transition Frequency logS 





Lower 


Upper 


Lvjnzj 




MT-T„ 
iNJ-l3 


z, 


J/Z, 


z 


z. 


J/Z, 


z 


9^ 799 '^Q1 ^90/1 n^i^ 


1 9/191 
— l.ZDzl 


IN Ms 


9 
z, 


J/Z, 




9 
z. 


7/9 
/ /z. 


J 


91 799^01 S7Q^^QA^ 

zj./zzjyio /yi^oOj 


1 Q01 


iM-13 


2, 


3/2, 


2 


2, 


3/2, 


2 


91 799fi191(\lC11^ 
Z J , / ZZO JZ jLW-^^ J J J 


^07Q 




2, 


5/2, 


3 


2, 


5/2, 


3 


23 722633644('33') 


0000 


NH3 


2, 


1/2, 


1 


2, 


1/2, 


1 


23.722634389(33) 


-0.7058 


NH3 


2, 


5/2, 


2 


2, 


3/2, 


2 


23.722673071(32) 


-1.2621 


NH3 


2, 


7/2, 


3 


2, 


5/2, 


3 


23.722675390(42) 


-1.3010 


NH3 


2, 


3/2, 


1 


2, 


1/2, 


1 


23.722679956(22) 


-1.0068 


NH3 


2, 


5/2, 


3 


2, 


5/2, 


2 


23.723887894(125) 


-2.0580 


NH3 


2, 


5/2, 


3 


2, 


3/2, 


2 


23.723928698(64) 


-0.9119 


NH3 


2, 


7/2, 


3 


2, 


5/2, 


2 


23.723929630(50) 


-0.7570 


NH3 


2, 


1/2, 


1 


2, 


3/2, 


2 


23.724686811(63) 


-1.0068 


NH3 


2, 


3/2, 


1 


2, 


5/2, 


2 


23.724691591(57) 


-0.7516 


NH3 


2, 


3/2, 


1 


2, 


3/2, 


2 


23.724732357(123) 


-1.7058 


i\i-l3 


-I 
J, 


Q/9 


J 


-I 
J), 


7/2, 


2 


91 8fi78n<;nQ4('99'> 


801 ^ 


IN 113 


J, 


7/9 


D 


J), 


5/2, 


2 


91 8^^781 /^^iQ1 /4'^^ 


— u.yooi 


lNJrl3 




J>/Z, 


■2 
J 


■2 


1/2, 


2 


91 8fi789zl'>1 1 9^ 
Z J.oD / oZM-J iO^^ i iZJ 




l>lrl3 


■2 

D, 


J/Z, 


■2 
D 


J), 


3/2, 


2 


91 8^789^808/90^ 


1 1 f^OO 

— ±.±ouu 


lNJrl3 


-I 
J, 


5/2, 


3 


J), 


7/2, 


4 


91 8fi84189zLl/A7^ 


1 OfifiQ 


IN 113 


T. 
J, 


7/2, 


3 


-J 


9/2, 


4 


91 868440144/94^ 


Q/19S 


iMl3 


J), 


3/2, 


3 


J, 


5/2, 


4 


91 8^^84zL^^9^^1 /71 


1 1 Q^S 


rMM3 


J, 


9/2, 


3 


■2 


11/2,4 


91 8A84'\01 96^^71 ^ 


— U.ozOo 


iNJrl3 


a 
J) 


3/2, 


2 


J) 


5/2, 


2 


91 870040070/91^ 




Mi-r 


'J 
J, 


3/2, 


3 


■2 
J), 


5/2, 


3 


91 8700ft487<;/4^;^ 
Zj.o /UUO^o / J^'I-OJ 


— U.O / 04 


iNJrl3 


^, 


7/2, 


3 


^, 


9/2, 


3 


91 8700^^7401 /'^1^ 




— U.oOO ( 


l>lrl3 


J, 


9/2, 


4 


3, 


11/2,4 


91 870070000/17^ 


Q1 ^0 

— u.y iou 


iNil3 


^, 


3/2, 


3 


^, 


3/2, 


3 


91 8701 97881 /90'\ 
Zj.o /UIZ / 00 H^ZUJ 


971 
— U.Z / ±0 


IN 113 


J, 


5/2, 


4 


J, 


5/2, 


4 


91 870190ft1ft/90\ 
Zj.o /UlZy010^^ZUj 


0000 

u.uuuu 


MT-T 
INJ-I3 




1/2, 


2 




1/2, 


2 


91 8701 10994/90^ 
Zj.o /Ul JuZZ'ti^ZUJ 


— U. / Ozi 


iNn3 


3, 


11/2, 4 


3, 


9/2, 


4 


93 870180310/60^ 


— n Q1 'M\ 


NH3 


3, 


9/2, 


3 


3, 


7/2, 


3 


23.870188326(66) 


-0.8557 


NH3 


3, 


5/2, 


2 


3, 


3/2, 


2 


23.870211423(62) 


-0.6059 


NH3 


3, 


11/2,4 


3, 


9/2, 


3 


23.871807418(23) 


-0.8258 


NH3 


3, 


5/2, 


4 


3, 


3/2, 


3 


23.871811331(20) 


-1.1938 


NH3 


3, 


9/2, 


4 


3, 


7/2, 


3 


23.871817154(22) 


-0.9428 


NH3 


3, 


7/2, 


4 


3, 


5/2, 


3 


23.871819253(22) 


-1.0669 


NH3 


3, 


3/2, 


2 


3, 


5/2, 


3 


23.872431263(30) 


-1.1500 


NH3 


3, 


1/2, 


2 


3, 


3/2, 


3 


23.872433410(22) 


-1.3541 


NH3 


3, 


5/2, 


2 


3, 


7/2, 


3 


23.872441474(22) 


-0.9651 


NH3 


3, 


7/2, 


2 


3, 


9/2, 


3 


23.872452975(39) 


-0.8013 
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Table S2: Absorption line parameters and l-cr uncertainties for the 8-component fiducial fit in 
Fig. 1. The optical depths, r, for the HCN and NH3 transitions are summed over all hyperfine 
components. For NH3, r is given in units of 0.01. 

Com- Redshift FWHM r r/lO"^ 

ponent [kms-^] HCQ+(l-2) HCN(l-2) NHgCl l) NH3(2,2) NHgCB 3) 

1 0.6845833(12) 2.80 ±0.40 0.0041^;^^^ 0.152l^:^i^ 0MT^^^^~0m^^^~0A9^^ 

2 0.6846071(12) 3.00 ± 0.20 0.4i;[|;}g OAbtojt O.OOlJ;? 0.37^^20 OM^oH 

3 0.6846205(34) 3.56 ±1.07 O.lStg'J^ 0.3ll[]:J? 0.61^^;^^ 0.32l°:f^ 0.23tiil 

4 0.6846473(133) 5.11 ±2.53 O.Sgt^J'as O-^SlS? 0.44l[j|j 0.57^0:3! 0.40tS 

5 0.6846568(9) 3.14 ±0.37 O-Hto'^o ^■'^'^-oil ^-^"^-oH ^-^^-oH 0-37;o:32 

6 0.6846730(19) 2.86 ± 0.98 0.59^111 O.eolo'sg l-26lo:64 l-14lo:65 0.85toi| 

7 0.6846879(38) 4.25 ±1.27 0.63lgJ^ 0.72lJ:^^ 0.86to:6o 0-32^0:30 O-l^lais 

8 0.6847024(138) 6.55 ± 2.01 0.25lg:i6 0-401^^4 0.43lg:^^ 0-46j:o:28 0-20lg:i4 
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